|cond-mat/9908398| , La Plata-Th 99/09, LPENS-Th-17/99, 

ETH-TH/99-20, SISSA 103/99/EP 



Doping-dependent magnetization plateaux in ]>merized Hubbard 

chains* 

D.C. Cabra 1 , A. De Martino 2 , A. Honecker 3 ^, P. Pujol 2 and P. Simon 4 

L Departamento de Fisica, Universidad Nacional de la Plata, C.C. 67, (1900) La Plata, Argentina. 
Facultad de Ingenieria, Universidad Nacional de Lomas de Zamora, Cno. de Cintura y Juan 

XXIII, (1832) Lomas de Zamora, Argentina. 
2 Laboratoire de Physique^ Groupe de Physique Theorique ENS Lyon, 46 Allee d'ltalie, 69364 

Lyon Cedex 07, France. 
^Institut fur Theoretische Physik, ETH-Hdnggerberg, 8093 Zurich, Switzerland. 
4 International School for Advanced Studies, Via Beirut 2-4, 34014 Trieste, Italy. 

(February 14, 2000) 

Abstract 



We study zero-temperature Hubbard chains with periodically modulated hopping at arbi- 
trary filling n and magnetization m. We show that the magnetization curves have plateaux 
at certain values of m which depend on the periodicity p and the filling. At commensurate 
filling n a charge gap opens and then magnetization plateaux correspond to fully gapped 
situations. However, plateaux also arise in the magnetization curves at fixed n between 
the commensurate values and then the plateau-value of m depends continuously on n and 
can thus also become irrational. In particular for the case of dimerized hopping (p = 2) 
and fixed doping we find that a plateau appears at m = 1 — n. In this case, there is still a 
gapless mode on the plateau leading to thermodynamic behavior which is different from 
a completely gapped situation. 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



Macroscopic quantum phenomena in strongly correlated electron systems in low dimen- 
sions are presently the subject of intense research. In particular, plateaux in magnetization 
curves of quantum magnets have recently received much attention, a central observation 
being that the plateaux occur at (typically simple) rational fractions of the saturation mag- 
netization. In one dimension, this is by now rather well theoretically understood in terms 
of a quantization condition that involves the volume of a translationally invariant unit cell 
(see e.g. [l|-|4j|). Various materials exhibiting plateaux have also been studied during the 
past few years, including e.g. a new candidate for a frustrated trimerized chain material 
||. One of the clearest examples is given by the low-temperature magnetization process of 
NH4CUCI3 || where one observes plateaux with 1/4 and 3/4 of the saturation magnetization 
thus demonstrating rationality even if it is theoretically still unclear why precisely those two 
numbers are observed in NH4CUCI3. 

The purpose of the present letter is to start a systematic investigation of the effect of 
doping on the quantization condition for the appearance of magnetization plateaux. We will 
show with examples that plateaux can appear in doped systems at magnetization values 
which depend continuously on the doping and are thus in general irrational. This obser- 
vation itself is not entirely new, but we believe our physical interpretation in terms of one 
class of excitations (say the up electrons) being pinned at commensurate filling while the 
filling of a different class of excitations (the down electrons) remains adjustable to be new. 
This interpretation also suggests that doping-dependent magnetization plateaux are prob- 
ably a generic phenomenon. One of the known cases is the one-dimensional Kondo lattice 
model where unpaired spins behave ferromagnetically thus giving rise to a spontaneous 
magnetization of a value controlled by doping. The other example is an integrable spin-S 1 
generalization of the t — J chain doped with (S — 1/2) carriers |J where, however, the 
appearance of plateaux is restricted to large magnetization values. 

As a first step towards a general understanding we study the effect of both a magnetic 
field and a periodic modulation of the hopping amplitude (p-merization) on a doped one- 
band Hubbard chain. One simple motivation is that for pure spin systems (corresponding to 
a half-filled Hubbard model) in a magnetic field, p-merized Heisenberg chains have turned 
out to be among the simplest examples M . In addition, structural modulations can give rise 
to dimerized coupling constants e.g. in chains and ladders ||, the organic (super) conductors 
|T0| (some of which come naturally at quarter filling) and the ferroelectric perovskites ]ll| . 

The one-band Hubbard model is not only used to describe realistic situations, but at 
least in one dimension it is also a useful model for technical reasons. For example, the 
Hubbard chain is exactly solvable by means of the Bethe Ansatz (BA) (T^j. This can be 
used as an input to a bosonization analysis [|l3H T5ll in order to study generalizations of this 
model (see e.g. P^|). 



In this letter we show that there are several different phases in the fi-h plane (see Fig. |l]for 
a schematic illustration of the cases p = 2 and p = 3; abbreviations refer to the corresponding 
regions of this figure): 
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i) If both quantization conditions^ 

|(n±m)eZ (1) 

are satisfied, both spin and charge excitations are gapful and hence all correlators 
decay exponentially at large distances (regions labeled by fixed n and m in Fig. [[]). 
In this case, we find plateaux in the magnetization curves since the presence of a spin 
gap is equivalent to the appearance of a plateau. 

ii) A charge gap ('CG') can open if the combination pn G Z of the conditions (|]) is 
satisfied. This includes the well-known charge gap at half filling (n = 1) and also the 
charge gap in the quarter-filled {n = 1/2) dimerized Hubbard chain {p = 2) |fl7f1 . 

iii) If only one of the conditions (|Ij) is fulfilled (solutions are indicated in the corresponding 
regions of Fig. |l|) and in addition doping n is kept fixed, a magnetization plateau opens, 
but one mode remains gapless (similar observations have been made in other systems 

19[). In contrast to the gapful magnetic behavior, charge transport remains 
metallic in this phase. 

The existence of these plateau phases is the main result of the present letter. A 
particularly appealing aspect of the plateaux predicted here is that they can appear at 
low magnetization (and thus at small magnetic fields) if the doping is chosen suitably. 

In the remaining cases, both spin and charge sectors are massless, leading to a Luttinger 
liquid ('LL'). 

II. HAMILTON! AN AND STRONG-COUPLING LIMIT 

To be concrete, we study the following model Hamiltonian 

L 

H =~Y< K X ) ( C Ul,aCx,a + cl a C x+ha ) 
x=l <J 

L L 

+U Y n x,\n X)i + /i (n x ,-[ + n X)[ ) 

x=l x=l 

h L 

~9 H K,t - "x,l) ( 2 ) 

Z x=l 

where t(x) = t if x 7^ mp and t(mp) = t' = t + 5. Here, and c are electron creation and 
annihilation operators, n Xi(T = c X(T c x ^ the number operator and a =f, [. \i is the chemical 
potential and h is a magnetic field. 

First we consider the limit t' = 0. Then the chain (0) decouples into clusters of p sites 
and one can use simple arguments in the spirit of : The number of up and down electrons 



18, 



n and m are normalized such that < n < 2 and Iml < 1. 
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on a p-site cluster must both be integer which is equivalent to imposing both conditions (HI). 
All these states are clearly fully gapped at t' = 0. Thus, they will remain fully gapped if 
one switches on a small perturbation t' > 0, only the transitions between these fully gapped 
states will soften. In fact, we will argue soon that these fully gapped states survive even 
until t' = t. 



III. BOSONIZATION 



Now we turn to a bosonization analysis starting with h = (for a related recent study 
of a dimerized Hubbard chain see |20|| ). For t' = t (5 = 0) and n ^ 1, the Hubbard chain 



can be represented by two bosonic fields with Hamiltonian O 



E | / dx [( W 2 + (Wl > ( 3 ) 

i=c,s 

where 4> = 4>r + 4>l and 9 = 4>r — 4>l- The spin and charge fields are given by (f> s = 771 (0t — 4>i) 
and (fi c = I (0| + (fii), where 0^ are compactified with periodicity — > 4>^,{ + V^" ^- The 
parameter £ and the Fermi velocities t> $ can be obtained exactly from the BA solution of the 
model in terms of U and \x. 

If we now turn on a non-zero but small 5, one can show using the continuum represen- 
tation of creation and annihilation operators in terms of the bosonic fields that the most 
relevant perturbations to the free Hamiltonian are given by 



Opert = A sin[A; + /2 + pk + x — \/7r£0 c ] x cos[v27T0 s ] + A' cos(A; + + 2pk + x — v47r£0 c ) , (4) 

where A and A' are proportional to 8 and k + = kp^ + kp,\, = nn. We can now study the 
values of n for which these operators are commensurate. In this way, if pn £ Z the A' term 
appears as a relevant perturbation, opening a charge gap. If we restrict further to pn/2 £ Z 
the A term is present as well and opens also a spin gap. 

We consider now the more complicated case of non-zero magnetization. The Hubbard 
chain in the presence of a magnetic field is also integrable |I3| , |T4|| . The large-scale behavior 



of the system is given by a two-field Hamiltonian like (|3D, but now the effective charge and 
spin fields are given by |T5[ 



! / 7 7—7 

' / Zj ss ss ^ cs 



det Z \ Z sc Z sc — Z cc 



(5) 

where Z is the so-called charge matrix given in whose entries can be obtained exactly 
from the BA solution in terms of U, /1 and h. These quantities determine the scaling 
dimensions of the vertex operators. 

For small 5 the perturbing operators read 

O pert = A sin[A; + /2 + pk + x - ^fix (Z cc (p c - Z cs (f) s )} 
x cos[A;_/2 + pk_x 

-\pn{{Z cc - 2Z SC )4> C - (Z cs - 2Z ss )(j) s )] 
+ A' cos(/c+ + 2pk + x - V47i(Z cc (j) c - Z cs (j) s )) , (6) 
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where now A;_ = kp^ — kp^ = irm. If all the operators are commensurate both degrees 
of freedom are massive since the perturbing operators are relevant. We thus have a mag- 
netization plateau with a charge gap. This is achieved when the two conditions (|Xj) are 
simultaneously satisfied, which generalizes the result for p-merized spin chains Q ton^ 1. 

When only one of these conditions is satisfied, say p(n + m) jl 6 Z, the Hamiltonian 
can be written as 



H 



+ 



dx 



+ A sin 2y/n(f)-i 



(7) 



Terms mixing derivatives of the up and down fields, which come from the U interaction, 
can be shown to be irrelevant using a treatment similar to Ref. [ 
the massive field 0-r and obtain an effective large scale Hamiltonian for 



One can integrate out 
with an effective 

Fermi velocity and Luttinger parameter K. This field is apparently massless but constrained 
to be in a particular topological sector if one imposes fixed filling (which seems to be natural 
from the experimental point of view). In a system of finite size L, this constraint takes the 
form Ln = -^(0t + 4>l)\o an d allowing only small energy fluctuations locks the quantity 



| q to be also constant 



see 



|21|j for more details). 



If this constraint were not imposed, Q would vary and the susceptibility would then have 
the standard form \ = K/(2tvv). On the other hand, if we impose Qo = L{n — m)/2, we 
obtain X = jz ((Q 2 ) ~ (Q) 2 ) = 0' indicating the presence of a magnetization plateau. Since 
the local part of the down sector in (|7|) remains massless, local correlation functions of the 
spin-down field decay algebraically and one can also show that the low-temperature specific 
heat is linear in the temperature. This situation is interesting since it combines a massless 
behavior with a magnetization plateau at p (n + m) /2 G Z for arbitrary values of n. 



IV. SMALL ON-SITE REPULSION 

A complementary derivation of these doping-dependent magnetization plateaux (at fixed 
n) is given by the following discussion of the limit of small interaction U. 

For the non-interacting case {7 = 0, the Hamiltonian (|2|) can be diagonalized easily and 
is found to have p bands e x (k). This is illustrated in Fig. ^| for the case p = 3 (the band 
structure for p = 2 is shown e.g. in Fig. 2 of |I7[). In the presence of a magnetic field h 0, 
the up and down electrons are subject to different chemical potentials (/i — h/2 and fi + h/2, 
respectively). A doping-dependent magnetization plateau then corresponds to a situation 
where one chemical potential (say for the up electrons) lies in one of the p — 1 band gaps 
while the other (for the down electrons) is in the middle of a band. Imposing the constraint 
of fixed filling n then requires a finite change in magnetic field in order to move the chemical 
potential for the up electrons into one of the bands next to the gap, leading to a plateau in 
the magnetization curve. However, the filling of the down electrons remains adjustable and 
one obtains a doping-dependent value of the magnetization on the plateau. 

A finite on-site repulsion U > leads to corrections to this non-interacting picture 
which can be treated to first order in U using standard quantum mechanical perturbation 
theory. We are interested in a situation where one band for say the up electrons is either 
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completely filled or empty. It is therefore sufficient to look at diagonal matrix elements of 
Hi = UY,x=i n x,t n x,i in the basis diagonalizing (|J) with U = 0, i.e. scattering processes (as 
would be present for half-filled bands) are absent in the present situation. This then leads 
to a simple shift of the critical fields h± for a doping dependent magnetization plateau. In 
particular, the plateau is non-vanishing also in the presence of interaction at least if U is 
small enough. 

For p = 2 this computation can be carried out explicitly (details will be presented in 



2jJ). One finds that the matrix elements of Hi are of a mean-field type, i.e. simple products 
of the filling of up and down electrons. The final result for the lower and upper boundaries 
of the m = 1 — n plateau for p = 2 is given by 



h± = ±\t - t'\ - ^ft 2 + t' 2 + 2W cos {{2n - 1) tt) + (n - l)U + 0{U 2 ) . (8) 

This leads to a plateau width h + — h- = 2\t — 1'\ + 0{U 2 ) which is independent of U at first 
order in U. The latter is a direct consequence of the mean-field form of the matrix elements 
of Hi and is specific to p = 2. 

We have also performed Lanczos diagonalizations of the Hamiltonian (|2]) for p = 2 and 
p = 3 at U = 3t, t' — 0.7t to confirm the overall picture. Here we just mention that the 
schematic Figs. [I] are in fact based on these numerical results, but otherwise postpone a 



detailed presentation to [21 



V. CONCLUSION 

To conclude, we have shown that there are two types of situations with plateaux in the 
magnetization curve of a p-merized Hubbard chain at arbitrary filling n: One where both 
pn and pm are integers and the complete excitation spectrum is gapful and a second one 
where only one of the conditions (0) is satisfied and part of the excitations remain gapless. 
If one works at fixed filling n, the latter leads to doping-dependent magnetization plateaux. 

Although we have studied p-merized Hubbard chains, doping-dependent magnetization 
plateaux should exist in more general situations and indeed, similar phenomena have been 
theoretically observed in other systems @,@]- 

From the experimental point of view, the most direct check would be a magnetization 
experiment on a p-merized chain with controlled doping. As we have mentioned earlier, 
materials realizing dimerized and trimerized chains do exist (also at non-trivial fillings). 
However, doping can be controlled better in a different class of quasi one-dimensional ma- 
terials, namely the high-T c related realizations of iV-leg ladders 0. In these materials, a 
doping dependence could be particularly intriguing since it could be used to push mag- 
netization plateaux into experimentally accessible field regions despite the large coupling 
constants. We therefore regard the magnetization process of doped iV-leg Hubbard ladders 
as an interesting field for further research. 
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